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Mechanisms of nonradiative recombination of electron-hole complexes in Cd(Mn)Se/Zn(Mn)Se 
quantum dots accompanied by interconfigurational excitations of Mn 2+ ions are analyzed within 
the framework of single electron model of deep 5rf-levels in semiconductors. In addition to the 
mechanisms caused by Coulomb and exchange interactions, which are related because of the Pauli 
principle, another mechanism due to sp-d mixing is considered. It is shown that the Coulomb 
mechanism reduces to long-range dipole-dipole energy transfer from photoexcited quantum dots to 
Mn 2+ ions. The recombination due to the Coulomb mechanism is allowed for any states of Mn 2+ 
ions and e-h complexes. In contrast, short-range exchange and sp — d recombinations are subject 
to spin selection rules, which are the result of strong Ih-hh splitting of hole states in quantum 
dots. Estimates show that efficiency of the sp-d mechanism can considerably exceed that of the 
Coulomb mechanism. The phonon-assisted recombination and processes involving upper excited 
states of Mn 2+ ions are studied. The increase in PL intensity of an ensemble of quantum dots in a 
magnetic field perpendicular to the sample growth plane observed earlier is analyzed as a possible 
manifestation of the spin-dependent recombination. 
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I. INTRODUCTION 

A great deal of attention attracted to quantum dot 
(QD) structures is due to the possible use of their quan- 
tum states in various fields of spintronics and for gen- 
eration and detection of light ji Semimagnetic (diluted 
magnetic semiconductor (DMS)) II- VI quantum dots are 
promising objects for these purposes because a high de- 
gree of spin polarization of electrons and holes can be 
achieved in relatively weak magnetic fields. 

However, incorporation of Mn ions into CdSe/ZnSc 
QDs substantially reduces the quantum yield of radia- 
tion as soon as the optical transition energy exceeds the 
energy of internal Mn transition ^2.15 eV. It is found 
that Cd(Mn)Se/Zn(Mn)Se QD photoluminescence(PL) 
is completely quenched due to nonradiative recombina- 
tion at relatively small Mn content x ~ 3 — 5%i 2 -~— 

The PL quenching is due to nonradiative recombina- 
tion of QD e-h complexes accompanied by excitation of 
Mn 2+ ions. The Coulomb interaction between 3d- and 
band electrons is usually considered as a reason for pro- 
cesses of energy transfer from e-h complexes to Mn 2+ 
ionsii^ 

The interest in possible mechanisms of energy trans- 
fer arose after observation of a strong increase in PL 
intensity of an ensemble of Cd(Mn)Se/Zn(Mn)Se QDs 
in magnetic field B || Oz reported by many authors. 2-6 
Simultaneous increase in the QD PL life-time was also 
reported£&i°i 

The explanation of the increase is based on suppres- 
sion of the spin-dependent nonradiative recombination 
of e-h pairs by magnetic field. According to the model 
proposed in Ref. 11 nonradiative exciton recombination is 



possible because of the direct exchange interaction be- 
tween band carriers and 3d electrons. This model pro- 
vides selection rules for the process: S' B = Sb, where 
Sb and S' B are Mn 2+ spin projections in the direction 
of magnetic field B in the initial and final Mn 2+ states 
respectively. Selection rules predict that the nonradia- 
tive recombination of bright excitons is forbidden for the 
states of Mn 2+ ions with Sb = ±5/2 whereas it is al- 
lowed for Mn 2+ ions with other Sb- Since the number of 
Mn 2+ ions with Sb = —5/2 increases with magnetic field, 
nonradiative recombination is suppressed. The rules cor- 
rectly explained the strong increase in PL intensity in 
Cd(Mn)Se/Zn(Mn)Se QDs but failed to explain its de- 
pendence of the direction of the magnetic fields To ex- 
plain the observed dependence a modification of selection 
rules was suggested in Ref. 2 S z + s ex , z = S' z , where S z 
is the projection of Mn 2+ spin in Oz axis, which is per- 
pendicular to the sample growth plane, s eXyZ — s e z + s h z 
is the projection of exciton spin (instead of the total 
angular momentum J z = j z + s|) in Oz. This sug- 
gestion is based on PL studies of Cd(Mn)Se/Zn(Mn)Sc 
and CdSe/ZnSe/ZnMnSe QDs, which revealed that hh 
and Ih hole states are strongly split due to strain and 
dimensional quantization! 12 i 13 The carriers in upper hh 
band are characterized by both the moment projection 
j z = ±3/2 and the spin projection s\ — ±1/2. The 
Coulomb and exchange mechanisms are assumed to obey 
the same selection rules ^ 

The analysis presented below shows, however, that 
spins of atomic configurations were not properly regarded 
in Ref.— and some correction of previous results is re- 
quired. Besides, the sp-d mixing specific to DMS was 
omitted while its contribution to the nonradiative re- 
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combination can be essential. The consideration based 
on the single-electron model of the 3d- level in a semi- 
conductor reveals that mechanisms based on the direct 
exchange interaction and sp-d mixing are subject to the 
spin selection rules In contrast, the Coulomb mechanism 
does not lead to the spin-dependent recombination. Es- 
timates show that efficiency of the sp-d mechanism can 
noticeably exceed efficiency of the Coulomb process in 
Cd(Mn)Se/Zn(Mn)Se QDs. The sp-d mechanism can ei- 
ther contribute to impact excitation of Mn +2 ions in- 
corporated into a II- VI semiconductor matrix, which is 
closely related to the nonradiative recombination. 

The paper is organized as follows. In the next section 
a model Hamiltonian of a DMS QD structure is consid- 
ered. Mechanisms of the nonradiative recombination due 
to the Coulomb interaction and sp-d mixing are analyzed 
in Sec. III. Phonon-assisted processes and energy transfer 
into upper Mn 2+ excited states are discussed in Sec. IV. 
The link between the model in question and the impact 
excitation of Mn 2+ ions is also discussed in Sec. IV. Spin- 
dependent selection rules for nonradiative recombination 
and their manifestations in experiments are considered 
in Sec.V. 



II. MODEL HAMILTONIAN 

The Hamiltonian of the structure containing one 
QD and one Mn 2+ ion has the form of an Anderson 
Hamiltonian)^— 

Hq = Hqd + Hmu + Hhyb + H Xl (1) 

where Hqd is the Hamiltonian of the QD electron sys- 
tem, Hhyb is the hybridization Hamiltonian, Hmu is the 
Hamiltonian of the Mn 2+ 5c?-shell and the Hamiltonian 
of the potential exchange interaction is H x . 
The QD Hamiltonian 

Hqd = H c + H v + H cv (2) 

contains Hamiltonians of the conduction H c = 
£ ^c+ Sz c^ and valence H v = V„, <o K. 
electrons. Here fi and v enumerate electron and hole 
states in QD. Conduction and valence electron states are 
KJ = c+ J°) and K*> = K 3 M, respectively. Elec- 
tron, hole, exciton, and trion states are eigen-states of 
Hqd- Hole states are related to valence electron states 
via a time- reversal operation = \<fZ*-j,)- 

The term H cv contains e-h Coulomb and exchange 
interactions. The Coulomb term mixes single-electron 
states and leads to the reduction of the e-h energy 
whereas the e-h exchange term splits e-h states with 
J = 1,2. Without loss of generality we retain only the 
exchange term and neglect the Coulomb interaction be- 
cause the confinement of e-h states in QDs is determined 
by the dimensional quantization and because its contri- 
bution to energy transfer is negligible as it will be clear 
in Sec.III. 



The dots under study have a very anisotropic lens- 
like form, i.e. the diameter D is several times larger 
than the height L. Photoluminescence studies of 
Cd(Mn)Se/Zn(Mn)Se and CdSe/ZnSe/ZnMnSe QDs re- 
veal that most of QDs have symmetry C 2v . These stud- 
ies show that hh and Ih hole states are strongly split due 
to strain and dimensional quantization! 12 ' 13 The e-h ex- 
change interaction splits e-h states with J z = ±1 ("bright 
excitons") and J z — ±2 ("dark excitons"). The gap be- 
tween bright and dark states is around 2-3 meVf^ 2 - The 
wave-functions of bright exciton in those structures are 
lV4c) — l/\/2(|l) ± | — I)), where | ± I) are exciton states 
with J z = ±1 and \tp^ x ) — | ± 2) are dark exciton states. 

The effect of the periodic lattice potential on band elec- 
trons is taken into account within the effective-mass ap- 
proximation, which may be used since we are interested 
in the properties of states near the bottom of the conduc- 
tion and top of the valence bands. In the effective mass 
approximation conduction electron states are ip e s (r) = 
F es Sz {r)S Sz +F ex Sz (r)X Sl +F ey Sz (r)Y Sz +F ez a ,{r)Z tm , 
where envelopes F es Sz (r),F ex Sz (r),F ey Sz (r),F ez Sz (r) 
are solutions of the effective mass equations. Here S, 
X, Y, Z are zone center Bloch functions of appropri- 
ate symmetries. The electron states on the bottom 
of the first band of dimensional quantization within 
symmetric quantum well of width L w can be found 
elsewhere. 17 Bloch amplitudes of free electrons with k || 
Oz that account for k 2 terms within the Kane's model are 
*4i /2 (fc) = (l-(hkp) 2 /6m 2 )S ±1/2 ±phk/(3m )[(l/E g - 
l/(E g + A))(X±iY) Tl/2 + (2/E g + l/(E g + A))Z ±1/2 }. ls 
Here p — (S\p x \X) — iPmo/h, P is the Kane parameter, 
mo is the free electron mass. 19 These results show that 
the admixing of Z functions to ip e s (k) is much larger than 
those of X and Y functions. 

The admixture of Ih to hh states takes place in QDs 
of C 2v symmetry and lower^ However, it is small and 
unimportant for our aims. Therefore, we assume that 
^ ±3/2 (r) = <; 3/2 (r) = ±l/V2F hhu (r)(X ± iY) ±1/2 , 
where basis functions from Refi^ 9 - are used. Thus, QD 
exciton states are | ± 1) = I ^±3/2) I ¥^1/2) and I ^ 2) = 

I ^±3/2>K±l/ 2 >- 

The Hamiltonian of the Mn 2+ ion in the Hubbard form 

is 

H-Mn = ^ fl mSz {e d + Ueffhm-sJ, (3) 
ms z 

where h mSr is the number operator of 5d-electrons with 
the orbital momentum to and spin projection s z , U e ff 
is the electrostatic repulsion energy between 5d-electrons 
occupying the same state to. We assume that Hmu acts 
only on d 4 , d 5 , d 6 configurations of Mn 2+ ions containing 
n=4,5, and 6 3d-electrons, respectively! 14 ' 20 The crystal 
field splitting of 5<i-shell is omitted in Hmu- When it 
is taken into account the Hamiltonian can be expressed 
as H Mn = J2i \A l )e l (A i \, where \A. k ) and e i7 respectively, 
are states and energies of d 4 -d 6 configurations. They can 
be found by means of the crystal field theory^ 
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The hybridization Hamiltonian 

Hhyb = H sc i + H p d (4) 
consists of Hpd = Em^.-jj^ ms.^j.t.k + ^- c -) 

and # sd = E rasrfs; (^ ms^^^C^' + ft - c -) terms - 

The sp-d hybridization results from the combined influ- 
ence of the potential of crystal ions and band electrons.— 
It is considered on the symmetry grounds as^ a result of 
action of an effective crystal field potential U cr , so that 
hybridization coefficients V sd mSzfJ . s > z = {d mSz \U cr \^s' z ) 
and V pd = (d m sjU cr \vj z ). In the case of substi- 

tutional Mn 2+ ions the crystal field mixes valence band 
states at k ~ only with ^-functions, which belong to 
t2 representation of the Td site symmetry^ - — The op- 
erator U cr also mixes 3d- and conduction band states at 

The direct s-d exchange term is 

H x = -£j(R A )S A ctsAs^^, (5) 

A vs z ,v' s' z 

where J(Ra) is the exchange constant, S> is the spin 
of the A-th Mn 2+ ion located at R\, § is the electron 
spin operator. The direct exchange is associated with 
the Coulomb interaction between band and 5c?-electrons. 
The direct p-d exchange is considered to be zeroi 14 i 15 ' 18 

Reduction of the PL quantum efficiency in semicon- 
ductor structures due to presence of impurities is usually 
related to the nonradiative recombination of e-h pairs be- 
cause of the Coulomb or exchange interactions between 
band and impurity states.— i 19 i 25 The Coulomb interac- 
tion is missed in the Hamiltonian in Eq.Q. In order to 
account for it the approach developed in Ref^ is used. 

The Anderson Hamiltonian in Eq. (p} , which was intro- 
duced phenomenologically into the DMS theory , 14 ' 15 can 
be obtained as a generalization of the Slater-Koster equa- 
tion for single-electron states in a crystal containing an 
impurity level (Ho + Hi mp — E)<j) e {v) = 0, where Hq is the 
single-electron Hamiltonian in an ideal crystal structure 
in the absence of the impurity, and Hi mp is the poten- 
tial of impurity 2 - The single-electron states are expressed 
via eigen-states 4>i{r) of H </? e (r) = F^^), where 
Fi are appropriate coefficients and summation is made 
over quantum numbers characterizing crystal states. 

The single-electron potential of perfect crystal struc- 
ture may be written as V° = £a «"( r ~ r a) + K°( r )r^ 
where V^°(r) is the potential due to band electrons, 
v A (r — R A ) is the potential of core electrons and nucleus 
at R A . In the presence of impurity the potential may be 
written as V = J2\ v *( r ~ r a) + V e (r), where Ra are 
new atomic sites. Therefore the perturbation is 

H imp = V - V° = £ A Mr - Ra) - «£(r - R A )) + 
U c (r) ~ 

^(, A (r-R A )-<(r-R A ))-£ ) p (6) 



x(R-R° A ) + L/ c (r), 

where U c (r) = V e (r) - V°(r). 

By neglecting lattice relaxation (Ra = R A ) the poten- 
tial of a substitutional impurity at ao may be written as 
H imp (r) = v 3d {r - a ) - v h (v - a ) + U c (r)^. whereas 
the potential of an interstitial impurity is H' imp (r — a^) = 
«3d(r — Sio) + U c (r). Here u/j(r — a ) is the potential of the 
host ion. Only substitutional impurities are discussed be- 
low since obtained results can be easily extended to the 
latter case. 

The Hamiltonian of the crystal structure with sub- 
stitutional impurity in the second-quantized form can 
be expressed as H = Hi + H2, where Hi = Hq + 
Hi m p is the single-electron Hamiltonian in the ba- 
sis {d mSz ,4>f, Sz ,4>J l j z }, that includes localized 5d-states 
d mSz . The term H2 contains interactions between 
single-electron states. The only interaction included 
in the Anderson model is the Coulomb repulsion be- 
tween electrons with the same m but opposite s z Hi = 
u effJ2 mSz ^ns z n m - Sz . Besides H QDl H hyb and H d = 

12ms e dd mSz d m s z the Hamiltonian Hi contains the scat- 
tering term 

H S cat — ^ t (U VSz ,u' s z c vs z c u' s z -\-U ll j z ,^> j z b ^^b^i j z ) 

H^H' ,v^v' ,s z ,j z 

(7) 

z,Mz C ts z b^j z + h.c.) 

usually omitted in the canonical Anderson model. Here 

U vs z ( Mz ).v's' z {^'j z ) = (vS z {Mz)\H im p\v' ' s' g (ji' ' j z )) . 

The Bloch wave- functions ^> VSz e lUr , where ^ VSz can 
be either X Sz , Y Sz , Z Sz or S s ,, are not orthogonal to 
dms z at k ^ 0. As a result anticommutators [d/x^c^] — 

(dms z \f% z ) # and [d mSz bt jz ] = (d ms Jflj z ) ^ be- 
cause slowly varying envelopes F eStex ,ey,ez,hh(^) can al- 
ways be expressed as £ k A(k)e lUr . We assume, how- 
ever, that overlap integrals are negligible because we are 
interested in band states near k = 0. 

To account for the sp-d Coulomb interaction between 
band and five 5d-electrons the terms originating from 
2i=i e2 /l r — r i| should be added to H^. The poten- 
tial U c (r) tends to compensate changes in electron den- 
sity caused by presence of the impurity and leads to 
the screening of the interaction between band and 3d- 
electrons. Within the framework of linear response the- 
ory H imp (q) = l/e e //(q)(ud(q)-«/i(q)), where H imp (q), 
Ud(q), Uft(q) are Fourier transforms of fl" imp (r), Ud{r — 
ao), it/i(r — ao), respectively. 26 Similarly, the screened 
Coulomb interaction between band and i-th 3d-electron 
is 

47re 2 f e lt i( r - r i) „ 

v sc {r - Ti) = -—. t / — r^2 d Q 8 

(27r) 3 J e e// (q)q 2 
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This expression is appropriate for the Coulomb pro- 
cesses of nonradiative recombination that are character- 
ized by small transferred quasi-momentum qu ~ 1/-D, 
q z ~ 1/L. Here D and L are characteristic lengths of 
QDs in the lateral and Oz-direction, respectively. Quasi- 
momentums within this range participate in the forma- 
tion of e-h states, which justifies the use of the lin- 
ear response approach.— Since the transferred energy 
Eq ~2.15 eV satisfies condition E ~ E g ^> Hujlo, 
the inert ionic system of the crystal does not contribute 
to the dielectric function so that the effective dielectric 
function is determined by band electrons. The processes 
caused by the long-range Coulomb interaction are char- 
acterized by q ~ so that for the Coulomb processes 

Contrary to the Coulomb processes exchange ones are 
characterized by a small interaction radius so that con- 
tribution of the processes with large transferred quasi- 
momentums can be noticeable. The explicit form of the 
screened Coulomb interaction is not required for the fu- 
ture analysis of the exchange mechanism so that the ex- 
pression W;j C (ri — r 2 ) is used in this case. 

The Coulomb interaction between band and 3d- 
electrons V sc = J2i=i w sc(r — r,) generates many terms of 
the form Vij^^af a^ai, where af(ai) can be either of 
rf ms iS ( d m S J ) c+ S2! (c /1Sa ),&+ 2 (6 I/ j ;a ). The terms containing 
only 3d- or band electron operators, i.e. those that do 
not couple 3d- and band states, are not considered here 
since the effect of these terms are already accounted for 
in H Mn and H QD . 



The other terms appeared because of 



V sc are 



analyzed by means of the approach developed in Ref. 27 
The analysis shows that some terms originated from V sc 
are already included in Hq. 

At distances |r — ao| as, where is the atomic 
Bohr radius, v sc (r — ri) ~ v sc (r — a ), so that the term 



ms z ,p,q 



*p,ms z 



(9) 



to 



compensate J2 P , q U p , q a+a q 



- a Q )\a q }a+a q that can be expressed as 



tends 

£^tp,q ^ 

J2p,q( a p\ z M n e 2 /e eff \r - a a \\a q )a+a q , where Z Mn is the 
net core charge that includes Mn 2+ ion nuclear charge 
and the charge of completely filled atomic shells. The 
expectation value (Y^ ms f l mSz — Zmu) is zer0 f° r an 
isoelectronic substitutional impurity^ The short-range 
difference between these potentials contributes to the 
potential scattering H scat and energies of the 3d and 
band states. 
Similarly terms 

^ (Kns s ,m' s' z .q.m' s' z ^m' s z d mSz d q ~\~ h.C.) (10) 

ras z ,m' s' z , q 



and 



^ ] {Ums z ,vs z d^ ns ^c lJSz + U mSzifl j z d^ nSz b fi j z + h.c.) 

(11) 



as well as 



ms z ,q,p 



(12) 



are included in Hhyb- The term in Eq. (fT2"]) tends to com- 
pensate the effect of crystal ions potential on the sp-d 
mixing.— 

Similarly the term 



E 



Vp,ms z ,p,m' s' z ^p^msz ^m' s' z 



(13) 



ms z ,m' s' z ,p 



tends to compensate the contribution of 



U, 



ms z .ms' ^ms z ,m' s' z ^ m s 



> s > d+ s d m < s i and the crystal ions 



£ 

core potential to energies of 3d and band states and the 
potential scattering of ^-states. 

By taking into account the Coulomb interaction the 
Hamiltonian in Eq.([T|) is transformed to 



H = Hq + H s 



U , 



(14) 



where Uq contains terms generated by V sc + Hi mp that 
are not included in Hq + H scat . They are responsible for 
various processes involving one, two or three 5rf-electrons 
and include those responsible for the nonradiative recom- 
bination. 

The Hamiltonian in Eq.(|14j) can be rewritten in the 
familiar form i 20 ' 29 



H = H, 



QD 



H 



Mn 



Hi. 



(15) 



where Hi nt is the interaction between the QD elec- 
tron system and the Mn 2+ ion. The hybridization term 
Hhyb in Eq. (fTfj) can be replaced by an effective scattering 
Hamiltonian i 22 ! 28 so that H can be rewritten as 



H = H, 



QD 



H 



Mn 



H hyb \I)(I\H, 



hyb 



Ei - Ei 



(16) 

where J5{ and Ej are energies of the initial \i) and 
intermediate \I) states, respectively. 

In addition to the terms H c e / / and H v e f / describing 
scattering of conduction and valence electrons on the 3d- 
shell the effective Hamiltonian contains the term H m i X 
responsible for the mixed scattering. 

As it is shown in Re L 22 i 28 the Hamiltonian 



H v 



I 



Hpd\I)(I\H p( i 
Ei — Ei 



(17) 



contains terms responsible for the potential scattering of 
valence electrons, contribution into energy of the 5c?-level 
and terms of the p-d kinetic exchange that can be ex- 
pressed in the Heisenberg form similar to H x in Eq.([5]). 22 
Analogous terms are contained in H c e ff. 
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The introduction of effective Hamiltonian H e f f is jus- 
tified in the ionic limit of the Anderson models It takes 
place when the broadening of the single electron 3d- 
states due to their hybridization with delocalizcd band 
states (QD states do not contribute to the broadening) is 
smaller than the energy gaps between the ground d 5 and 
excited d A , d 6 configurations.™ In this limit the ground 
and excited states of Mn 2+ ions are characterized by cer- 
tain spins S = 5/2 or 3/2, which is obviously the case 
here. 

The initial state, which we denote as [d eh], contains 
the Mn 2+ ion d 5 in the ground states and the e — h 
pair. Its energy is Ei = E (5) — E v + E c , where E (n) 
is the energy of d n configuration of the Mn 2+ ion, E v 
and E c are energies of electron states at the bottom of 
conduction and the top of the valence bands. Interme- 
diate states |J) are of two kinds: |/ + ) corresponding to 
the virtual absorption and |/ _ ) corresponding to the vir- 
tual emission of one electron. The dominant contribu- 
tion to energy transfer comes from the virtual levels with 
minimal energy gaps Ei - E ± . The energies of those 
states [d 6 h] and [d 4 e 2 h] are E+ = E Q (6) + U eff - E v , 
and E~ = E (4) + 2E C — E v , respectively. The values of 
energy gaps Ei — E^ are found in Appendix B. 

By using commutational relations and the fact that 
virtual levels involve only d e and d 4 configurations the 
term 



H,, 



E 



H pd \I)(I\H sd + H sd \I)(I\H pd 



Ei - Ei 



(18) 



energy migration between Mn ions. The influence of pairs 
and triads on energy transfer is considered in Sec V. The 
discussion of the dynamics of Mn-Mn energy transfer is 
out of the scope of current manuscript. 

Thus, TJq + H m i X can be considered as a perturbation 
of the Hamiltonian 



H' = H, 



QD 



H 



Mn 



He 



(21) 



may be expressed a; 



,28 



where H ex is the sp-d kinetic exchange term. The po- 
tential scattering is dropped as it does not directly con- 
tribute to the nonradiative recombination. 

The Coulomb and direct exchange interactions leading 
to the nonradiative recombination are contained in Uq. 
In the following section we consider contributions of Uq 
and H m i x to the recombination separately and estimate 
relative efficiencies of the Coulomb and sp-d mechanisms. 



III. MECHANISMS OF NONRADIATIVE 
ENERGY TRANSFER TO 3D-SHELL 

A. Nonradiative recombination due to Coulomb 
interaction 

The nonradiative recombination due to the Coulomb 
interaction is usually considered as a transition of the 
first order i 7 ' 8 ' 25 The only term in Uq leading to such a 
transition is 

5 

V t = ^( V rn,n, v , m >dt ls ^C+ s ,d m , s b v jz + 



E 



(19) 



V , d + r + b d , 



h.C. 



(22) 



E 



K 



^ + h.c), 



where coefficients are 



"ms.m's^W; — v sd ms^us z V pd m 's' z uj z ^ _ 



(20) 



and K mix = K+ Szm , <v ^ + # m ^ m , s , The second 
term in Eq. ([T9l leads to the interband potential scatter- 
ing whereas the first one contributes to the nonradiative 
recombination as it is shown in the following section. 

The Hamiltonian given by Eq. (|16[) can be generalized 
to the case of many non-interacting Mn ions which cor- 
responds to the limit of small Mn content when the in- 
teraction between Mn ions can be omitted. The antifcr- 
romagnetic Mn-Mn coupling appears within the Ander- 
son model in the fourth order of perturbation series and 
becomes important starting from x ~ 0.05. The Mn-Mn 
interaction leads to the formation of Mn-Mn pairs, triads 
and clusters. The antiferromagnetic Mn-Mn interaction 
leading to the cluster formation affects the dynamics of 



where V, 



V n . 



m,fi,m' ,v 



(d mSz ,fis' z \v sc \d m , Sz ,isj z ) and 



m' s z I ■ 



Contrary to the potential scattering where the single- 
electron term Hi mp tends to compensate the effect of 
sp-d Coulomb potential, Hi mp does not participate in 
the nonradiative recombination because matrix elements 
between initial and final states including different 3d- 
configurations vanishes. 

The rate of resonant energy transfer from the pho- 
toexcited QD into Mn 2+ M-shell is given by the Fermi's 
golden rule: 



1 



2tt 1 



Y^\M fi \ 2 5{E t ~E f ), 



(23) 



if 



where N$ is the number of initial states^ 

The initial state of the system |i(5/2, S z , G)) = 6 
Ai(5/2, S z )i/j e x{G)\0) consists of 5 5d-electrons and the e- 
/jpair while the final state \ f(3/2,S' z )) = 4 fi(3/2,^)|0) 
is the excited state of the Mn 2+ ion. Here ip ex (G) = 
4>ex(G)\0) is the exciton state, where G denotes a set 
of quantum numbers characterizing it. In the absence 
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of magnetic field G is a representation of the symmetry 
group of exciton state. For instance, it is E for bright 
and A X ,A 2 for dark excitons in quantum dots and wells 
of Did symmetry. The exciton creation operator ip ex (G) 
is expressed via c^ Sz and b V j z operators. In strong mag- 
netic field B || Oz exciton states are characterized by 
J z . Operators 4 fi(3/2,S' 2 ) and 6 ii(5/2,S" z ) are sums 
of products of five d+ s ^ operators. Wave-functions of 
the many-electron 4 Ti(3/2, S* z ) = 4 T\(3/2, S Z )\Q) and 
6 Ai(5/2,S z ) = 6 Ai(5/2, S z )\0) states in various approx- 
imations of the crystal field theory can be found in 
literatures^ 

The matrix element of nonradiative recombination 
Mfis^AG = (f(3/2,S' M )\V t \i(5/2,S,,G)) reduces to a 
sum of Coulomb and exchange integrals. The Coulomb 
integrals lead to dipole-dipole energy transfer from QDs 
to Mn 2+ ions as it is shown in Appendix A. The anticom- 
mutation between 3d and band electrons does not affect 
Coulomb matrix elements so that energy transfer due to 
the Coulomb interaction can be understood as a result of 
emission and absorption of virtual photons by QD states 
and Mn 2+ ions. The crystal media can screen this pro- 
cess but cannot interfere with it.— The dipole transition 
between 4 Ti(3/2) and 6 Ai(5/2) states is spin-forbidden 
and, therefore, admixing of Mn 2+ excited states with spin 
3/2 to 6 Ai(5/2) is required. 

Robbins and Dean considered the dipole-dipole mech- 
anism as the dominant in nonradiative exciton decay ac- 
companied by intraionic excitations. 7 Similar point of 
view is expressed in Re fi 25 ' 30 ! 32 where excitation of Er 
ions in Si accompanied by optically forbidden inter- 
configurational transitions of Er is considered. How- 
ever, efficiency of the exchange mechanism can substan- 
tially exceed that of the Coulomb one when the latter 
is spin-forbidden, which is usual in atomic and molecu- 
lar systems i 33 ^ 34 Analyzing the process of deexcitation of 
Mn 2+ ions in the presence of free carries, in other words, 
the reverse process in respect to the impact excitation 
(Auger process), Allen concluded that the Coulomb in- 
teraction underestimated it as much as by two orders of 
magnitude and, therefore, the main mechanism of energy 
transfer was related to the exchange interaction rather 
than the Coulomb one.— The idea about domination of 
exchange mechanism in excitation of Mn 2+ ions is widely 
accepted now i 11 ' 35 ' 36 The exchange mechanism of energy 
transfer in atomic and molecular systems is subject to 
the Wigner spin conservation rule, which states that it 
is allowed if the total spin of the interacting system is 
conserved^ 

The spin conservation rule was used by Nawrocki et al. 
to derive selection rules for the exchange mechanism. 11 
However, in order to apply it the authors of Ref. 11 simpli- 
fied the studied system so that they neglected any spin- 
orbit coupling in the system and consider the transition 
\ 6 A 1 {5/2)(p e s ^) | 4 T!(3/2)^, ) instead of exciton recom- 
bination. Exciton states in QDs are not completely char- 
acterized by spin due to the strong spin-orbit interaction 



in the valence band and e-h exchange interaction. How- 
ever, spin selection rules for the exchange mechanism can 
be derived from the analysis of exciton and Mn 2+ spin 
functions. 

The matrix element of the nonradiative recombination 
due to Coulomb interaction in the first-quantized form is 
Mji = (f\V sc \i), where \i) = ^A^h^S^)), \f) = 
| A( 4 Ti(3/25 z )V'o), V'ex is the exciton many-electron func- 
tion, Vo is the ground state of the crystal, and A is the 
antisymmetrization operator. 

The exciton states are sums of products of spacial 
and spin functions corresponding to spins s ex = 0,1 and 
spin projections s ex z = for bright and s ex z = ±1 
for dark excitons because the wave-functions of 
bright exciton ground state can be expressed as 
^(ri.ra) = F hh0 (r x )F &0 (r 2 )/V2(-(X + ^y) 1 S 2 a ^1 /3 2 ± 

{X-tY)iS 2 0ia 2 ) = -F hlM {r 1 )F e0 {r 2 )/V2{X 1 S 2 {a 1 p 2 T 
fi\ot 2 ) + iY\S 2 {a\[3 2 ± f}\a 2 )) whereas the 
dark exciton wave-functions are ip^ x (ri,r 2 ) = 
-F hh0 (n)F e0 (r 2 )/V2{X + iY) 1 S 2 a x a 2 , and 
V4(ri,r 2 ) = F hm {r x )F e0 {r 2 )/y/2{X + iY) x S 2 ^ 2 . 
Here a and /3 are s z — 1/2 and s z = —1/2 spin states, 
respectively. 

Unlike the Coulomb processes the antisymmetrization 
between 3d- and band electrons is crucial for the ex- 
change mechanism. The wave-function of the final state 
|/) is that of the crystal in the ground state ipo contain- 
ing one Mn 2+ ion in the excited state 4 T!(3/2, S z ). The 
final state is a sum of Slater determinants constructed 
from N + 5 single-electron wave- functions where N is the 
number of valence and conduction electrons of the crystal 
including two Mn 4s electrons that become delocalized. 
Each Slater determinant can be presented as a sum of 
products of spatial and spin functions, where spin func- 
tions are basis functions of the irreducible representation 
of the permutation group of N + 5 electron spins corre- 
sponding to certain squared spin S 2 and spin projection 
S^, whereas spatial functions are basis functions of the 
representation conjugate to the spin representation^ 7 - It 
means that the final state certainly contains spin func- 
tions with S = 3/2. 

The initial state \i) obviously contains Slater determi- 
nants with spin projections S z = ±5/2 and S z = ±7/2, 
which means that corresponding spin functions have 
S=5/2 and 7/2. They can also have 5/2 - 1 = 3/2 con- 
tained in 5/2 ± s ex . It means that the nonradiative tran- 
sition is allowed for bright and dark excitons. Similarly 
it is allowed for X~ trions. When the exchange mech- 
anism is allowed it necessarily leads to the conservation 
of the spin-projection S z + s ex z because of properties of 
the exchange integral V mtiUm > . 

The exchange matrix elements do not allow 
parametrization and hardly be estimated in a sim- 
ple manner. They can be small or large because the s-d 
exchange constant is about a quarter of the magnitude 
of the p-d kinetic exchange constant^ whereas contri- 
bution of the direct exchange into the p-d exchange is 
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usually assumed to be zeroi 14 ' 18 

The use of the Fermi golden rule is correct if the matrix 
elements of perturbations H mix and V t are much smaller 
than energy difference between any Mn 2+ configurations 
and e-h QD states. It is valid for the sp-d and Coulomb 
mechanisms as it is shown in Appendices. It is natural to 
expect that this is also valid for the exchange mechanism. 

The hybridization Hhyb makes it possible Coulomb 
processes via virtual states involving d 4 and d 6 config- 
urations because of the second-order term 



E 



(i\U \I)(I\H hyb \f) (i\H hvb \I)(I\Uo\f) 



E t - Ei 



Ei-Ei 



(24) 



Its contribution to energy transfer is small first of 
all because matrix elements Hhyb if /(Ei — E{) and 
Hhyb n/{Ei — Ei) are much smaller than unity as it is 
shown in Appendix B. Besides, the Coulomb matrix ele- 
ments between \i) \f) and intermediate levels |/) involv- 
ing d 4 and d 6 configurations with spin S — 2 is spin- 
forbidden. 



B. Nonradiative recombination due to sp-d mixing 

The sp-d mixing is responsible for the nonradiative re- 
combination that can be understood as a result of suc- 
cessive hopping of the electron(hole) and hole (electron) 
in the 5c?-shell. The hopping of a valence electron in and 
out of the 3d-she\l is the reason for the kinetic p-d ex- 
change interaction. The hopping of conduction electrons 
becomes possible because of the admixture of valence 
band states to conduction ones as k ^ 0. It is respon- 
sible for the reduction of the s-d exchange interaction 
reported in Refj 1 ^ The sp-d mechanism was proposed by 
Schmitt-Rink et at as a mechanism of excitation of rare- 
earth ions incorporated into a semiconductor matrix* 38 
It is related to the term 



Hmix /J (E-mixd mSz d nl f s' z ^ v j z C^ls z ~\~ h.C.)^ 



(25) 

The p-d coupling coefficient for the valence elec- 
tron ground state, v = 0, is V pd mSs0 j 3 =±3/2 = 
(d m sjUcr\Vo ±3 / 2 ) = TFhho{ao)/V2(d ms jU cr \(X ± 
iY)±i/2), where a is the position of Mn 2+ ion. The 
s-d coupling coefficient for the conduction electron 
ground state fi — is V sd ms'^Oa* = (d ms >jU cr \ipQ S J = 
F ex o( a o)(dms'jU cr \X Sz ) + F ey o(a.o)(d ms 'JU cr \Y Sz ) + 
F ez0 (a )(d ms i z \U cr \Z Sz ) since (d ms > z \U cr \S Sz ) = because 
of the symmetryi 14 ' 15 i 22 ' 23 

Due to strong anisotropy of QDs under study the main 
contribution to V s d is because of quantization along Oz. 
The coefficient of s-d mixing for the electron ground state, 
v = 0, in a QD can be expressed as: 



V sd 



ms*', ..0s z 



F eza (aL )(d m \U cr \Z)S s 



(26) 



7nz u s' .s z 

at k close to the center of Brillouin zone. Coefficients 
Vpd ms z 0j z and V sd ms z o Sz are related to the p-d hopping 
amplitude V p ° d as it is shown in Appendix B. 

The matrix clement of recombination of the bright ex- 
citon l 6 ^!^)^ = ±1) | 4 ri(^)> calculated in Ap- 
pendix B is: 



M„ 



16 



a(5 z )— F, iM (ao)F ez0 (a )|y; d 



|2 



(27) 



U, 



eff + e d 



+ 



E r E, 



Ed 
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where a(S z ) are certain coefficients found in the Ap- 
pendix B. The factor ^7/2 appears because of admixing 
of valence band states to conduction ones. It is impor- 
tant that such a transition conserves the Mn 2+ ion spin 
projection S z = S' z . 

The p-d exchange constant appearing in the hh ex- 
change Hamiltonian H ex = — /3/3jS(ao)<5(r — ao) is given 
by the expression : 14 ' 15 ! 22 



32 1 
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(28) 



This expression is obtained neglecting the crystal field 
splitting of Mn 2+ states i 14 ' 20 This result is reproduced 
with due regard for crystal field splitting of the initial 
stated It is natural to expect that neglecting the crys- 
tal field splitting in matrix element M m i X in Eq. (|27p also 
does not noticeably affect the result. The physical mean- 
ing of sp-d mechanism is illustrated in Fig. 1. The pro- 
cess of kinetic p-d exchange interaction is shown in Fig. 2 
for comparison. 

The sp-d mechanism of nonradiative recombination 
can be much more effective than that due to the Coulomb 
interaction. Estimates presented in Appendices A and B 
provide the ratio of recombination rates of the Coulomb 
and sp-d mechanisms for CdMnSe/ZnSe QDs averaged 
over distribution of Mn ions Rcoui/Rmix — 10~ 2 . 
The envelope F es (r) reaches maximum at the center 
of QDs, whereas F ez (r) reaches maximum at the QDs 
boundaries,— therefore the ratio Rcoui/ Rmix depends on 
the distribution of Mn ions within the e-h pair localiza- 
tion volume. The ratio Rcoui/ Rmix is expected to be 
larger than Rcoui/ Rmix in CdMnSe/ZnSe and smaller in 
CdSe/ZnMnSe QDs. 

Although the sp-d mechanism was proposed for exci- 
tation of rare-earth ions in a semiconductor— it better 
matches our case. For instance, ^/-orbitals of Er ions in 
Si are located much below the top of valence band (~ 10 
eV),— which results in weak p-f and s-f hybridizations. 
Besides 4f-orbitals are closely located to the ion core so 
that their interaction with bands electrons is weak. In 
contrast, the 5rf-level lays at 3.4 eV below the top of va- 
lence ban d 14 ' 15 and resonantly mixes with valence band 
states. 
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The dipole-dipole mechanism is an analogue of the 
Forster mechanism of energy transfer between atoms or 
molecules whereas the exchange mediated mechanism is 
similar to the Dexter one^ The mechanism related to 
the sp-d mixing can also be associated with the Dexter 
mechanism, where H m i X plays a role of the effective ex- 
change interaction. 

The contribution of the direct exchange mechanism to 
the nonradiative recombination can hardly be estimated 
in a simple manner so that relative contributions of M ex 
and M m i x is unknown. However, both mechanisms lead 
to spin-dependent energy transfer. Possible manifesta- 
tion of such a process is discussed in Sec. V. 

It is worth noting that the sp-d mechanism can con- 
tribute to the impact excitation of Mn 2+ ions. The ex- 
citation via either intra- or interband electronic transi- 
tions is possible . 25 ' 38 The impact excitation by means of 
the intraband transition becomes allowed when the en- 
ergy of an optically or electrically excited electron with 
respect to the bottom of conduction band exceeds 2.15 
eV. The impact excitation can be understood as a result 
of capture of such an excited conduction electron on the 
3d-shel\ because of the s-d mixing with successive escape 
of one of 5<i-electrons with the opposite spin also into the 
conduction band. Ayling and Allen state, however, that 
the efficiency of excitation of Mn 2+ via interband transi- 
tions is much larger than that via intraband ones.— ^ 



IV. PHONON-ASSISTED PROCESSES AND 
ENERGY TRANSFER INTO UPPER MN 2+ 
EXCITED STATES 

The resonant processes discussed above assume equal- 
ity of the e-h pair transition energy and the energy of in- 
traionic excitation Eq ~ 2.15 eV. Experiments show that 
the photoluminescence quenches even when transition en- 
ergy of the ground state of e-h pairs in QDs substantially 
exceeds Eq&&- There are two mechanisms, which allow 
such processes: the phonon-assisted recombination and 
the nonradiative transition in upper Mn 2+ states. 

The model of phonon-assisted energy transfer from 
a semiconductor crystal to rare-earth impurities is pro- 
posed by Yassievich et ali& It is based on the single con- 
figurational coordinate scheme describing the electron- 
phonon interaction of impurity electrons with a phonon 
mode of the crystal. According to the model nonreso- 
nant recombination is possible due to emission of mul- 
tiple phonons. The probability of the phonon-assisted 
transition is given by: 

R nr = f ^ EE \M fi \ 2 I(N ph )6(E g -E -N ph hu ph ), 

1 if N ph 

(29) 

where N p h = (E g — Eo)/fkJ p h- At low temper- 
atures ,fcflT <C hujph, the coefficient I(N p h) ~ 
(SH) Nph e~ Nph /Nphl, where Sh is the Huang-Rhys fac- 
tor, which determines the strength of electron-phonon 



coupling^ 9 . According to the recipe of Refi^ 5 - summation 
over Nph is replaced by integration and leads to the ex- 
pression: 

W/*) = - = ^f^Ei^i 2 ( 3 °) 

T nr h TUjJph N *-rf 

if 

The probability Rf^(Kfi) has the maximum at Kfi — 
Sh at large coupling strength Sh > 1 and quickly de- 
creases with Kfi. The upper limit of the efficiency of the 
phonon-assisted recombination can be estimated by the 
summation over all Kfi. 

The single configurational coordinate picture is often 
used for the analysis of temperature dependence of PL 
lines of Mn 2+ internal transitions and the Stocks shift be- 
tween Mn 2+ absorbtion and emission spectra . 24 ' 35 ' 40 The 
structures under study are characterized by more than 
200 meV wide Mn 2+ PL line 2 - that can be attributed to 
strong interaction of the 3d-shell with crystal vibrations 
(Sh ^ 1). Assuming that the energy of the vibrational 
mode is of order Iiujlo one obtains from T nr ~ 10~ u 
s for a Cdo.85Mno.15Se/ZnSe QDs with reasonable values 
of D ~ 35A and L ~ 12A . 

Although this value of r nr is not far from measured 
PL decay times t = 20 - 80 ps of CdMnSe/ZnSe QDs 
with close parameters at B — T^l small intensities 
of the longitudinal optical phonon (LO-phonon) repli- 
cas of Mn 2+ PL lines in DMS materials indicate weak 
coupling between the 5d-shell and LO-phonons^ 2 - Sim- 
ilarly these structures are characterized by a weak vi- 
bronic interactional In this case recombination processes 
into upper excited Mn 2+ states can dominate because 
nonradiative energy transfer is possible not only into the 
lower excited state 4 Ti but also into other excited states 
4 T 2 , 4 Ei, 4 Ai located 0.2-0.7 eV above 4 T X . The large in- 
homogeneous broadening of Mn 2+ PL lines makes non- 
radiative transitions possible in a wide range of energies 
higher than 2.15 eV. The excited 4 Ti, 4 T 2 , 4 A l7 4 E states 
have spins S = 3/2, so that relaxation into 4 Ti state is 
expected to be very fast. 

V. SPIN-DEPENDENT SELECTION RULES 
AND INCREASE IN PL INTENSITY IN 
MAGNETIC FIELD. 

The strong increase in PL intensity of an ensemble of 
Cd(Mn)Se/Zn(Mn)Se QDs in the magnetic field B || Oz 
is explained as a result of suppression of the exciton 
nonradiative recombination because of depopulation of 
A 1 {5/2,S Z ) states with S z > -5/24*11 The results of 
Sec. Ill reveal that the Coulomb process is allowed for 
S z = ±5/2 and it does not contribute to the increase in 
PL intensity. 

The exchange and sp-d mechanisms lead to conser- 
vation of the total spin projection S z + s ex ^ z = const. 
Recombination of bright excitons (s ex , z = 0) requires 
AS^ = 0, whereas the recombination of dark excitons 
(s e x,z = ±1) is possible when AS Z = =pl. 
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The selection rules predict that the nonradiative re- 
combination of bright excitons involving S z — 5/2 is 
spin-forbidden, whereas it is allowed for the dark exci- 
ton with J z = 2(s ex . z — 1) and forbidden for that with 
J z = —2(s ex ^ z = —1). The strong increase in PL inten- 
sity requires slow relaxation of bright excitons to J z = 2 
dark states. Phonon-assisted recombination and energy 
transfer into upper excited Mn 2+ levels do not break the 
rules. 

However, the selection rules are not valid for e-h com- 
plexes formed by hole states with large admixture of Ih- 
states, which is the case for bulk DMS materials of cu- 
bic symmetry. In contrast, the selection rules should be 
valid for highly anisotropic crystals of wurtzite structure 
such as CdMnS, which is in agreement with the results 
of photoreflectance measurements reported by Nawrocki 
et alM 

Recent studies of the optically detected mag- 
netic resonance give additional arguments in support 
of spin-dependent energy transfer from ensemble of 
CdMnSe/ZnSe QDs to Mn 2+ ions and its dependence 
on the direction of magnetic field4£ 

The observed increase in QD PL intensity indicates 
that spin-dependent sp-d and exchange mechanisms dom- 
inate over the Coulomb one. This assumption correctly 
explains magnetic field and temperature dependence of 
the PL intensity of QD ensemble. The following empirical 
expression is found to fit the experimentally observed in- 
crease in PL intensity with magnetic field B || Oz: 1(B) — 
A/(l + Cp(B/T)), where p(B/T) ~ a + (3 e -*E(B)/k B T j 
a and j3, A and C are constants. 2 '- This well-known ex- 
pression describes the temperature dependence of the in- 
terband PL intensity in the presence of nonradiative re- 
combination centers. The magnetic field dependent ac- 
tivation energy is AE(B) = A e ^ + Ajw n , where A e ^ = 
Ms(3e + 9h)B is the energy of Zeeman splitting of e-h 
states. The parameter a is the probability of nonradia- 
tive recombination independent of magnetic field while /3 
is the probability of nonradiative recombination involv- 
ing the S z = —3/2 level^iii Results of Ref^ indicate 
that the ratio 1(B)/ 1(0) ~ (1 + /3/a) can reach 10 2 in 
high magnetic fields AE(B) /ksT 3> 1, which is possible 
if P <C a and the relaxation of the bright excitons into 
lower dark states is slow. 

Since the radiative life-time of 4 Ti(3/2) -> 6 A 1 (5/2) 
transition is substantially longer than the time of non- 
radiative recombination, fast saturation of 1(B) is ex- 
pected because only a few tens of Mn 2+ ions can inter- 
act with the localized e-h pair. The derivation of 1(B) 
implies, however, that the reservoir of Mn 2+ ions is infi- 
nite. This fact can be related to the fast energy diffusion 
within ensemble of Mn 2+ ions42 The dependence 1(B) 
very well describes experimental results^i 4 ^ which sup- 
ports the assumption. 

The QDs considered in Ref.— ^ are assumed to be neu- 
tral meanwhile considerable amount of QDs in the en- 
semble can be negatively charged due to a n-type back- 
ground doping almost always present in II- VI materials. 



The selection rules predict that X~ j z — ±3/2 states 
are involved in nonradiative recombination at B = 0. 
In magnetic field B || Oz, which polarizes both trion 
and Mn 2+ states only j z = +3/2 trion participates in 
the nonradiative recombination whereas j z — —3/2 trion 
does not. 

At Mn content of x = 1 — 2% the lower X~ trion 
state in magnetic field B j| Oz is that with j z = 3/2. 13 
Negatively charged QDs do not lead to the increase in 
PL intensity when the relaxation of the j z = —3/2 to the 
j z = 3/2 state is fast. 

However, measurements of X" trions in individual 
CdSe/ZnSe/ZnMnSe QDs reported in RefJ^ show sur- 
prising increase in the intensity of upper a~ component 
of the trion PL in magnetic field B || Oz. This result 
can be understood on the ground of the proposed model 
if the relaxation j z = —3/2 to j z = 3/2 is slow. The 
two level model, which takes into account the nonradia- 
tive recombination, relates the time of nonradiative re- 
combination t,„. to the time of spin relaxation r s from 
the upper trion state: r nr (B) < r s /2(l - e - A M n /kT^ 
where Amu = ^B9MnB is the Zeeman splitting of near- 
est Mn levels. This is probably because the relaxation 
—3/2 — > 3/2 requires changes Aj z = 3 so that it can be 
slow in QDs under study because of strong Ih-hh splitting 
and the splitting of Mn 2+ spin states in magnetic field. 
Negatively charged QDs, therefore, can contribute to the 
increase in 1(B). 

In the Voigt geometry (B _L Oz) magnetic field aligns 
Mn 2+ and electron spins opposite to the direction the 
field whereas hh hole moments remain directed along Oz. 
In addition, all dark | ±2) and bright | ± 1) exciton states 
mix at B 1 Oz. According to the selection rules the 
nonradiative recombination of dark excitons J z = 2 is 
allowed at any B, and, therefore, no increase in PL in- 
tensity is expected in the Voigt geometry, which is in 
agreement with the results presented in Rcf. 2 The non- 
radiative recombination of trions in the Voigt geometry 
is also allowed because of mixing of ±3/2 states. 

The Mn-Mn interaction has been omitted in the fore- 
going discussion, meanwhile it causes formation of Mn 
clusters among which Mn-Mn pairs and triads are most 
important. 44 Number of the pairs quickly increases with 
increase in Mn content and reaches maximum at x of 
several percents. The pairs are characterized by the rel- 
atively large antiferromagnetic coupling energy, about 1 
meV for nearest neighbors in CdMnSe and ZnMnSc, 14 
so that they remain coupled even at high magnetic field 
of about 12 T. Such pairs do not participate in the Zee- 
man splitting of band states as the total spin of the pair 
is zero. In contrast, the pairs can participate in radia- 
tive and nonradiative recombinations as those processes 
involve excitation of individual Mn 2+ ions in the pairs. 
The Mn-Mn coupling energy is still too weak to modify 
or mix 3d 5 configurations of individual ions within the 
pair as the energy gap between 5d-configurations is sev- 
eral eV. Therefore, there are not many reasons to expect 
strong changes of either Rcoui or the time of intraionic 
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optical transition l/r r . There are no reason to expect 
strong changes of M mix and M ex due to the pairs and 
triads formation. 

Experimental results reported in Ref4^ revealed, how- 
ever, shorter PL life-time of d-d pairs line than that 
of individual Mn 2+ ions. It can be assumed that 
the pair formation modifies Mn 2+ ion environment and 
changes covalent coupling strength, which, in turn, de- 
creases T r . Consequently, increase of Rcoui and the ratio 
Rcoui I ' Rmix with Mn 2+ content can be expected as R m ix 
linearly depends on x. The Mn-Mn pairs do not lead to 
the increase in 1(B), because magnetic field does not 
align spins of coupled Mn 2+ ions. Contrary to the pairs, 
Mn-Mn-Mn triads are characterized by nonzero average 
spin and therefore can contribute to the increase in 1(B). 

The sp-d and exchange mechanisms are short-range 
ones, they are effective within the e-h pair localization 
volume. On the contrary, the dipole-dipole mechanism is 
a long-range one. It may be important in CdSe/ZnMnSc 
QDs and obscure spin-dependent effects. At distances 
between the CdSe/ZnMnSe QD edge and the Mn 2+ ion 
core Rd much larger than the QD dimensions Mcoui de- 
creases as R2 6 i 31 i 33 Thus, the ratio Rcoui/ Rmix increases 
with Rd and the dipole-dipole mechanism can dominate 
at large distances. Assuming that the energy transfer 
between Mn ions is fast one can expect larger increase 
in 1(B) for CdMnSe/ZnSe neutral QDs with x < 0.1 
and the smaller one for CdSe/ZnSe/ZnMnSe neutral QDs 
which is in agreement with experimental results^" 41 - At 
large x > 0.1 the fast saturation of I(B)/I(0) or even 
decrease can take place. This is not valid in general case 
because the ratio 1(B)/ 1(0) depends on various parame- 
ters such as distribution of Mn 2+ ions, the ratio of neutral 
and charged QDs and probably considerably vary from 
sample to sample. 

The considered mechanisms of nonradiative recombi- 
nation are effective not only in DMS QDs but also in 
bulk and other low-dimensional structures. DMS II- VI 
materials are characterized by close parameters U e f / and 
E v — ed but slightly different E g , V s d and V p d- With the 
use of parameters reported in Refi 14 i 16 one can find that 
the contribution of the sp-d mechanism to energy trans- 
fer in CdMnS and CdMnTe structures is comparable to 
that in CdMnSe ones. 

The proposed model deals with direct energy transfer 
from e-h complexes to Mn 2+ ions. In bulk DMS materials 
and low-dimensional structures other processes via var- 
ious intermediate states such as defects, impurities and 
surface states are possible. The mechanisms considered 
above are involved in energy transfer in those cases too. 



VI. CONCLUSIONS 

Mechanisms of the nonradiative recombination of e- 
h complexes in Cd(Mn)Se/Zn(Mn)Se QDs accompanied 
by the intraionic excitation of Mn 2+ ions are analyzed 
within the framework of the single-electron model of deep 



5d-level in semiconductors. Together with traditional 
mechanisms related to the Coulomb and exchange inter- 
actions between 3d- and band electrons another mecha- 
nism caused by the sp-d hybridization is considered. Esti- 
mates of matrix elements reveal that the efficiency of this 
mechanism can considerably exceed that of the Coulomb 
mechanism and its contribution to nonradiative recombi- 
nation can be significant or even dominant. Mechanisms 
of energy transfer from neutral and negatively charged 
QDs to Mn 2+ ions due to the sp-d mixing and direct ex- 
change interactions are subject to the spin selection rules 
S z + s e x z — const in magnetic fields B || Oz, whereas the 
Coulomb mechanism does not obey them. These rules 
are because of the strong hh-lh splitting of hole states in 
Cd(Mn)Se/Zn(Mn)Se QDs. It is shown that nonradia- 
tive recombination remains efficient even when the fun- 
damental energy gap E g substantially exceeds the energy 
of the lower Mn 2+ internal transition most probably be- 
cause of energy transfer into upper Mn 2+ excited states. 
The proposed model indicates that the increase in PL in- 
tensity depends on the Mn 2+ content, distribution of Mn 
ions, QD dimensions so that its magnitude can strongly 
vary from sample to sample. 
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Appendix A: Matrix element of nonradiative 
recombination due to Coulomb interaction 

The only term in Uq leading to the transition from 
the initial \i) = 6 A 1 (5/2, S z )tp ex (G)\0) to final |/> = 4 
2i(3/2, <S Z )|0) states in the first order is 

5 

Vt = ( V m,^m' ^d+ Sz C+ s ,K^d ml .^ (Al) 

flips' tfziSz m,m' 

where V m ^m' ,v = (4j, , ^s' z \v sc \d m > Sz , vj z ) and 

Vm,iL,i>.m' (d< mSz , [IS z \v sc \v j z , djn' Sz ) . 

The Coulomb mechanism of nonradiative recombina- 
tion originates from the first term in Vt- The screened 
Coulomb potential is 

e 2 47re 2 f e l ^~ r ^ , s 

^ sc (r 2 - ri) = ; 1 = -— T / d q 

eoo|r 2 -ri| e 00 ( 27r) 3 7 q 2 

(A2) 
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The dominant contribution to the matrix element of 
long-range Coulomb potential comes from the domain 
|r 2 — a | » |ri — a | ~ Ob. By using the multipole 
expansion of the Coulomb potential 



v sc (r 2 - ri) = 



eoo|(r 2 - ao) - (ri - a ) 
(ri - a )(r 2 - a ) 



(A3) 



r 2 - a 



the matrix element V, 



:oo|r 2 - a |° 

is parameterized as 



v; 



e 2 e 2 

too ^oo 



/ n ( r 2 - ao) i . , 



(A4) 



|r 2 - a | 

The Coulomb operator can be rewritten 

&S Vf Coul Emm'sj ^™^"i's z ^ms^m's! 

E^^V^j/Js',^ + h - c - and tne Coulomb 
matrix element (f\V t Coul\i) is equal to e 2 /eoo(/|dP|i), 

where d = Y. m , m ' A rns z m> sAsA™' is the °P" 
erator of the intraionic dipole transition, and 

The matrix element 

Prtuj. = (Kl ,. 8 ", 80 ?, !^) = (A5) 



r 2 - ao 



(2^ 



d 3 



q q 



r 



(^4 



= iq(r 2 -ao) 



The wave-functions of the ground state of conduction 
electron,/^ = 0, are (r|0s' 2 = ±1/2) = (r|c^ 1 y 2 |0) = 

V5±i /2 (r) = F es0 (r)u C!±1/2 (r) = £ k A=(k)e' ikl V ±1/2 (r), 
where u C)=F i/ 2 (r) = S Tl/2 - 

The weak admixture of valence band states to con- 
duction ones is unimportant in calculation of P us > V j z 
and can be omitted. The wave-functions of the 
ground state of valence electron, v = 0, arc (r\0j z — 
±3/2) = ^ ±3/2 (r) = <r|fe+ 3/2 |0) = F^(r)« 0)±3/9 (r) = 

E k 4„(k)e ikr M„ :± 3/ 2 (r), where u V} ± 3/2 (r) = T{X + 
iY)±i/ 2 /V2 are periodic parts of Bloch functions of 
conduction and valence electrons, respectively. Wave- 
functions of the bright exciton are | ± 1) = c 4 ^ ; 2 &±3/2|0)- 



0s' z 0j z 



-*™*Zo ( a o)Fhho (ao)R^i. because 



Pi o s 'Oj e = (0s' x \r 2 j\0jz) 



4ir 



d 3 q qiQj 



(A6) 



Here the relation (niUj) = 47r/3<5y for the angle inte- 
gration is used, where n is the unit length vector. The 
vector 

1 



R 



r u 



Q 



d 3 7 



, {r)ru v , jlt (r) 



S2o 



(A7) 



where is the unit cell volume, can be R±i/2±3/2 = 
^fl/V2(r x cv ±iv v cv ) =^fr cv l/y/2(e x ±ie y ). Reree x ,e y 
are Ox and Oy unit vectors, r x cv — (S\r\X), r y cv = 
(S\t\Y), r cv = (S\y\Y) = (S\x\X) = (S\z\Z). 

The Coulomb matrix element between the ini- 
tial \i{5/2,S z ,G)) = 6 i 1 (5/2,5,)^« B (G)|0) and 
the final \f(3/2,S' z )) = 4 Ti(3/2, S z )\0) states is 
Mcoui-fS^s* = -e 2 /eood r 7 s^s,(0|P|i/ ; e2;(G)), where 
d r 7 s:,s z = (Tij(S' z )\d\Ai(S z )) is the matrix element of 
the 6 Ai(S z ) -^ 4 Tx 7 (S' z ) dipole transition. The transi- 
tion is spin and parity forbidden so that the spin-orbit 
interaction and admixture of p states to d m s z should be 
taken into account otherwise d r ^s' ,s z = 0. According 
to the model of Boulanger et ai4^ the optical transition 
6 Ai — > 4 Ti 7 takes place by means of spin-orbit interac- 
tion via intermediate state 4 Ti(3/2) that allows for the 
p-d mixing. 

Because of the symmetry d rx ~ / siS z r x cv = 
dry -fS' z s z r y cv = d rz jS ^ Sz r z cv and (d r 7 s' z S z v cv) 2 = 
r 2 „|d r -fS'S z \ 2 ■ The matrix element of transition 
6 A 1 (5/2,"s z )i'ex(Jz = ±1)|0) -> 4 fi 7 (3/2,^)|0) is 

4ne 2 l^j 

Mcoui jS' z s z — ± -r e so( a o)-F)i/io(ao)d r 7 s^,g z r el , , 

too 2y2 

(A8) 

which means that the Coulomb process reduces to the 
dipole-dipole energy transfer. 

The above result can be reproduced 
if the matrix element Mccml -yS' ,S Z = 
( 4 T lj (3/2,S' z )\V sc \ e A 1 (5/2,S z )i Jex ) instead of 
(f\V t coui\i) is considered, where V sc (r,r 1 ...r 5 ) = 
Y]j—i v sc (ri — r). Indeed, the matrix element of transi- 

tion \*Ai(5/2,S M )il>M = ±1)) -> | 4 T l7 (3/2,^)) can 
be expressed as 



Coul ■ySi.Sz 



k,k' 



(A c (k)A* v (k')u: tjz (r) 



e i(k-k0r4 T * (3/2) S ' z )\V sc \Uc,s z {vfA, (5/2, S z )) = 



47re 2 1 



eoo (2^)3 ^ 



^A c (k)A:(k')e j(k - k ' )ao (A9) 



x / d3 <l^ J s,Kj,k'( ( l) b iS' z ,s z {<i), 

^ where 6^,5, (q) - E< 1 1? 7 (3/2, S^e-^C"-") 

6 Ai(5/2, S z )d 3 ri ~ — zqd r 7 s'S-, because 

EJ 4 Tl 7 (3/2,^)r I 6 ^i(5/2,5 z )d 3 r l 

dr 7 s^s z - The matrix element J Sj kj^k'(q) = 

(^^(r)!^^-" 1 ' 4 -^^-^)!^,^^)) = i(k-k'+q)yR* ajV 

With that, the matrix element Mcoui 7 S' s z becomes: 

Mcoul iS' z s z — ± ^ d r 7 5,s^ r c „A c (k) 

e °° k.k' 
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xe lkao e" 



lk ' a °A* v (k') 



2V2 



Appendix B: Matrix element of nonradiative 
recombination due to sp-d hybridization 



47re 2 1 ~F ^ 

± F es0 (ao)F hh0 (ao)d r 7 scs z r CT — =, (A10) 

£oo 2y2 

Matrix elements d r 7 g' $ z are related to the radiative 
life-time of 4 Ti — s> 6 A\ transition: 



1 S' Z S Z 



3A^ 4 c 3 1 _3 2 
' "4^^^ "2^' (A11) 



where r r is the radiative time, Af, = (2 x 5/2 + 1) is the 
degeneracy of the state 6 Ai(5/2, S z ), n r = 



t 3 is 

the optical refraction index, and Eq ~ 2.15 eV is the 
transition energy; (D r /as) 2 ~ 10 -4 for CdMnSe since 
r r = 200 /xs4^ The factor 3/2 appeares because only 
drx ■yS' z s z and d ry ^s'S x components are involved in the 
nonradiative recombination; the factor 1/3 is because of 
degeneracy of the Ti 7 (3/2, S z ) state. 

Summation over Mn 2+ ions interacting with the bright 
exciton gives the rate of the transition (\ 6 Aiip ex ) — > 

l 4 ri»: 



2n rA ^2 1 
Rcoui = Y— 

n Ni \ eooCiB 



a%D. 



8m c a 2 B E g 



(A12) 



!^ S o(aA)| 2 |iWaA)| 2 ^ - E f ), 



where as is the atomic Bohr radius, and A enumer- 
ates Mn 2+ ions. In derivation of this equation rela- 
tions r = ihp/m E g , p = iPm /h, and P 2 = 3E g (A + 
E g )h/ (2(2 A + 3E g )m c ) ~ E g h 2 /2m c are used, where 
m c = 0.13rno is the conduction electron mass in CdSe 
or ZnSe; 19 P is the Kane constant. Since A/E g < 0.3 in 
ZnSe and CdSe crystals we assume that A/E g =0. 

After averaging over positions of Mn 2+ ion, the rate 
becomes: 



— 2?r, x9 1 

Rcoul = -j-^-kYNmuq 



D r 

(IB 



2 al 



V e V h 



r) e r]h (A13) 



8m c a 2 B E g 



5(E l - E f ), 



where V e nA is the effective volume of electron (hole) local- 
ization and coefficients ??e(ft) = Jdms d 3r \F e ( n ) (r)\ char- 
acterize penetration of the electron(hole) wave-function 
into the DMS layer. 

The presented analysis generalizes the approach devel- 
oped in RefjL£i2&, where the interconfigurational optical 
transitions are considered as single-electron ones. 



The Hamiltonian responsible for the nonradiative re- 
combination via sp-d mixing is 



E±rn.ix, — 



where 



E 

m,7n' ,v,fi,s z ,s' z ,j z 



(Bl) 



Kn 



sd ms z fis 



v* 



1 



1 



Ei - E+ E, 



E 
(B2) 

and Ei is the energy of the initial state whereas E + 
and E" are energies of intermediate states. For calcu- 
lations of the matrix element we use wave-functions of 
the ground and excited states in the spherical approx- 
imation similar to the approach of Schriffer in the cal- 
culation of the kinetic exchange constant. 22 We choose 
the initial state of the Mn+ 2 ion to be |Ai(5/2, -5/2) = 
^m=-2^m -1/2 1^)- States with 
tions can be obtained from this 
the step-up S+ = J2 m Kn 1/2^™ -1/2 an d step-down 
S- = J2m d m -i/2 d m 1/2 spin operators: |5/2,5 Z ) = 

where C is the normalization 



other spin projec- 
state by means of 



CS 5 _ /2 - s *l£ 



-2 d+ m 1/2 |0>, 



m 1/21 

constant. Particularly, C = l/v5 for S x = ±3/2, and 
C = 1/VlO for S z = ±1/2. 

We assume here that Hund's rules are valid. Thus, 
the final state is G 4 (3/2) excited state of Mn 2+ ion 
13/2,-3/2) = 4 1/2 4_ 1/a nU-i< 1/a |0>. Creation 
and annihilation operators df n , d m satisfy standard 
commutational relations. Although ((p e \d m ) ^ and 
(<p v \dm) 0, Kane basis functions are orthogonal to 3d- 
states, so that the nonorthogonality is negligibly small. 

The initial state for the recombination of 
bright exciton |1) 



C - 1/2 ^3/2^^-^ 



Tf 2 d + 

±1 m=-2 u m 1/2 



ii /2 & 3/2|0) is \i) 
|0), whereas the final 
state is |/) = d+ 1/2 d+_ 1/2 C'S 3 / 2 - s 'U} n= _ 1 d+ 1/2 |0). 

Since V p d couples valence band states with 5d-states 
of t2 representation it is useful to express the former 
via states of tetrahedral symmetry: d±2 = \j\[2(d ev — 
d t2 f) oc (x±iy) 2 , d±\ = =pl/ \/2(dt 2 (±idt 2 ^) oc y(x±iz), 
do = d eu oc (3z 2 — r 2 ). Here £, £, rj are basis functions of 
the i 2 representation, u and v are basis functions of e rep- 
resentation of the tetrahedral group £21 Matrix ele ments 
of the s-d coupling are |V; d ±2 | = \V sd ±1 | = \V pd \y^y/2. 
Transitions involving the state do are forbidden because 
it does not mix with valence band states, i.e. V sc t q = 0. 

By means of the relation V p d = AV^/^/TTq, where 
V? d is the hopping amplitud o 14 ' 22 the matrix element 



pd 

of the transition \Ai(±3/2)J z = 1) 
(|4i(T3/2)J, = -1) -> |Tx(t3/2))) is 



1^1 (±3/2)) 



M„ 



16 



a(±S/2)—F ez o( ao )F hho ( ao )\V p ° d 



|2 



(B3) 
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X \U eff +e d -E c + E v -e d ) (2) ' 

where No is the number of unit cells per unit volume and 
a(±3/2) = —1/y/E. By means of spin-up operators and 
wave-functions of initial and final states the coefficient 
a(±l/2) = -y/3/10 for the transition \A 1 (±1/2)J Z = 
1) ->■ |Ti (±1/2)) can be found. 

By using ip e s (k) for free electrons the coefficient 7 
is evaluated as 7 ~ h 2 /2m c E g L 2 . Smaller value 7 ~ 
h 2 /2m c E g L 2 /(l + h 2 /2m c E g L 2 ) can be found for the 
symmetric quantum well of width L with infinitely high 
barriers by using functions i y 9^_ 1 y 2 (r) presented in Refill. 
Since Ek = 5e d - E v + E c , E+ = 6e d + U eff - E v , E~ = 
Ae d + E e , U eff = 7.0 eV then E t - E~ = -(E v - e d ) ~ 
-3.5 eV, Ei-E+ = -{U eff + e d - E c ) ~ -1.5 eV. Val- 
ues of U e ff , E v — e d and V® d ~ —0.6 eV are taken from 
RefJ£ 

Taking into account Eq. (|B3p the probability of nonra- 
diative recombination is as follows: 

Rmix = = j-r ^2 $Z -T~\M m i X \ \ 2 5{Ei — Ef) 

A tf 

(B4) 

Here A enumerates Mn 2+ ions. The recombination rate 
of the bright exciton averaged over the Mn 2+ ions dis- 
tribution can be obtained in the manner described in 
Appendix A. 



The ratio of rates, corresponding to the Coulomb and 
sp-d mixing processes is: 

Rcoul = 1_( e 2 V fa%\ 2 (DA* 
Rmix ZyV^aBj \n J \a B J 

\UV pd J \a B J 

where U = l/(e d + U eff - E c ) + 1/(E V - e d ), £oo ~ 6, Q 
is the volume of unit cell. 

The parameters of U e ff, E v — e d , and V® d are not in- 
dependent but related to each other via the exchange 
constant /3. These parameters are obtained by means 
of theoretical and numerical analysis of Zeeman splitting 
of band states and photoemission spectroscopy dataJ^ 
The experimental results of Zeeman splitting as well as 
photoemission data are reported with accuracy of two 
significant figures. The best correspondence between ex- 
perimental and theoretical results are found with values 
Ueff = 7.0 eV, E v - e d = 3.5 eV, V° d = 0.6 eV for 
CdMnSe i 14 ' 16 Taking into account uncertainty of the pa- 
rameter e d + U e ff — E c related to the fact that E g de- 
pends on the Mn content and strain distribution within 
the sample we conclude that the accuracy of our estimate 
is better than an order of magnitude. 
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FIG. 2: The scheme illustrates processes of the kinetic p-d 
exchange interaction between hh 1 3/2) state and Mn 2+ 3d- 
shell. Two different paths of transitions via virtual interme- 
diate "donor" d 4 and "acceptor" d 6 states are shown. 
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FIG. 1: The scheme illustrates the process of the nonradia- 
tive recombination of the dark exciton J z = 2 in a magnetic 
field B || Oz. Two different paths of the process via virtual 
states contributing to (|27p are shown. The sp-d mechanisms 
of exciton recombination can be described as the result of suc- 
cessive hopping of the electron(hole) and hole (electron) in the 
3d-sheR. 



